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ABSTRACT 

Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static 
and dynamic operations of semiconductor lasers. It is well known that the value of diffusion coefficients affects the 
threshold pumping current of a semiconductor laser. At the same time, the strengt o carrier 1 usion process 
is expected to affect the modulation bandwidth of an AC-modulated laser. It is important not only to investigate 
the combined DC and AC effects due to carrier diffusion, but also to separate the AC effects from that o e 
combined effects in order to provide design insights for high speed modulation. In this presentation, we apply a 
hydrodynamic model developed by the present authors recently from the semiconductor Bloch equations. The mode 
allows microscopic calculation of the lateral carrier diffusion coefficient, which is a nonlinear function of the carrier 
density and plasma temperature. We first studied combined AC and DC effects of lateral carrier diffusion by s u ying 
the bandwidth dependence on diffusion coefficient at a given DC current under small signal modulation. The results 
show an increase of modulation bandwidth with decrease in the diffusion coefficient. We simultaneously studied the 
effects of nonlinearity in the diffusion coefficient. To clearly identify how much of the bandwidth increase is a result o 
decrease in the threshold pumping current for smaller diffusion coefficient, thus an effective increase of DC pumping 
we study the bandwidth dependence on diffusion coefficient at a given relative pumping. A detailed comparison of 
the two cases will be presented. 
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1. INTRODUCTION 

It is known that inhomogeneity in semiconductor lasers influence their operations . 1 - 11 Such non-uniformity can 
result from both extrinsic and intrinsic sources. For example, interface roughness and disordered alloying profile cause 
additional scattering for carriers, and thus reduce their mobility and diffusivity. Additionally, current spreading or 
crowding due to inadequate design can lead to excessive thermal budget, mcrease m injection current threshold, an 
reduction in slope efficiency for the active device. In particular, non-uniform carrier distribution and the associated 
diffusion process in the lateral quantum well (QW) layer have been shown, both experimentally and numerically, 
to affect dynamic behaviors , 2 modulation response , 3 mode dynamics 4 and selection, beam quality and optoelectnc 
property , 6 threshold current , 7 etc. Using a density-dependent diffusion coefficient, Sarzala and Nakwaski and the 
present authors 11 concluded that an average constant value of the coefficient is sufficient for reliable simulation 
results, which will be shown in this work the need for quantitative modification. 

Despite all the efforts, no correlation has been established for the effects of lateral diffusion on the DC operation 
and AC-modulation performance. In other words, how much of the change in the AC-modulation performance 
because of lateral diffusion is derived from the corresponding change in the DC characteristics. The answer o 
this question deepens our understanding of the static and dynamic nature of the effects of lateral diffusion and 
improves our insights for device design. Statically, lateral diffusion increases the threshold current and reduces the 
slope efficiency of an active device under influence, which in turn alters the device’s relative pumping level. According 
to the rate-equation description with a linearized gain model and no diffusion , 13 the corresponding spectral response 
under AC-modulation scales with the relative injection level and thus changes in full correlation. So it becomes 
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an apparent question: Will this prediction change in a description with a fully nonlinear gain model in addition 
to inclusion of diffusive process? In this work, we present such a model, select the relatively simple single mode 
gam-guided VCSEL case, and show that within the limit of our study, it can be concluded that the effects of lateral 
diffusion on AC-modulation (termed AC effects later, in comparison with the DC effects) are purely of static nature, 
which indicates that the AC effects are a result of changes in the DC operation characteristics, i.e., the threshold 
injection current and the slope efficiency. Definitely, possible changes are rich and non-trivial for cases such as a 
laterally-coupled strip laser (refer to the special session for this specific case in this proceedings) or a VCSEL featuring 
transverse mode dynamics. Additionally it is unambiguously shown that the nonlinear diffusion coefficient can be 
substituted by an equivalent constant coefficient as far as both DC and AC-modulation performance are concerned. 
A quantitative criterion is established to locate the constant coefficient. 

The paper is organized as follows. In the following section, theoretical aspect of our work will be summarized. 
The complete set of equations for our physical model is given together with expressions for the nonlinear density 
diffusion coefficient. Then in Section 3 , we first present the microscopically calculated diffusion coefficient as a 
function of plasma density and temperature, followed by numerical results and discussions for both DC operation 
and AC-modulation. We intend to identify the static and dynamic nature of the lateral carrier diffusion effects 
on AC-modulation by first applying a scaling scheme and then lineshape analysis for the spectral response under 
AC-modulation. Finally, the paper is concluded with a summary and some concluding remarks. 

2. PHYSICAL MODEL AND EQUATIONS 

The present work is based on a hydrodynamic model we developed recently 10,14,15 to describe the intra-QW 
carrier dynamics of a semiconductor laser. Next, we summarize the relevant equations and briefly explain the 
physics involved. First of all, our model explicitly considers differentiation between carrier types, which then can 
be conveniently reduced to a single-component description — the ambipolar diffusion (AD) regime — under certain 
approximations. This is the regime considered in this work. Furthermore, we neglect plasma heating effects in this 
study m order to reduce the associated complications. Their effects are topics of future study. 

In the AD regime, Poisson equation is eliminated due to charge neutrality and the whole set of equations for 
the description of the plasma density N , the complex transverse laser near field envelop £, and the macroscopic 
polarization P, within the QW plane f, is listed below: 


d t N + dr ■ (-D NN dfN) = -7 nN + Hi - ^3(P*£) , 

4 h y 




e 

P - k£ , 


d t Py = {-Ti (N, Ti ) + i[u c - a> x (Ar,2))]}P - ieotiA^N ,Ti)£ , 
P = P 0 + Pi ; P 0 = e 0 e bX o(N, T t )£ , 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


where J is the injection current with a spatial pumping profile and quantum efficiency rj, L m is the aggregate width 
of the active region, and T is the optical mode confinement factor. Shorthand notation for the partial derivatives 
is used, such as d t N = dN/dt. Other notations in the above equations are standard, for example, u> c is the central 
frequency and v s is the group velocity of light in the background material. 

In the above set of equations, the macroscopic polarization and laser field equations are treated following Ref. 
16. The polarization P is decomposed into two parts, electronic and background, as expressed by Eq. 4, where the 
background contribution Po is given in a constitutive relation as well. The electronic part is governed dynamically by 
Eq. 3. The parameters— the effective background susceptibility *o, the gain bandwidth T : , the peak gain frequency 
u>1 ’ 311 ^ Lorentzian oscillator strength Ay — are fitted to microscopically computed results as a function of carrier 

density N and lattice temperature I}. 16 The laser field is obtained after integrating over the assumed longitudinal 
mode profile. 


Distinctively in our model, the ambipolar diffusion coefficient P ; vA r is obtained microscopically by taking into 
account many-bodv effects within the Hartree-Fock approximation where carrier and polarization scatterings are 
treated under second-order Born approximation. The coefficient is found to be contributed by both free carrier and 



many-body parts. Specifically, it is dependent upon the scattering rates of electrons and holes with the longitudinal 
optical (LO) phonons in a semiconductor QW as listed below: 

D nn = (d N W + Nd N 5e g ) / (m e7 Lo + m h ^ 0 ) , ( 5 ) 

where W is the plasma thermal energy, 6c g is the bandgap renormalization owing to the Coulomb interaction between 
electrons and holes, m e , h and 7 are the masses and scattering rates of electrons and holes. 


3. SIMULATION RESULTS AND DISCUSSIONS 

Now we present and discuss our numerical results based on the above equations, which have taken into account 
the nonlinearity in both the material gain and the diffusion coefficient in the AD regime. The simulations axe run for 
a single mode gain-guided MQW GaAs/InGaAs VCSEL operating at an <n c of 980 run. The injection current profile 
is uniformly circular with a diameter of 7.5 /xm in the middle and smoothly reduced to zero at 8.25 pm fo owing 
a squared cosinoidal waist shape. The cavity length L is 144 nm and L m is 4x9 run -36 nm which yields a T 
of 0.25. The other VCSEL parameters and simulation details can be found m Ref. 17. Except the carrier lie ime 
( T = l/nrjy) of 0.25 nanosecond, all rates are computed microscopically and fitted as a function of l^asma ensi y 
N and lattice temperature Tj for our simulations. Furthermore, we use an injection current vit in 0 a o\e 
threshold current. Numerical results in this work are for near field at the center of the device, which is representative 
for single mode operation. The lattice temperature T h as well as the plasma temperature T p , is assumed to be a 
constant of 295 K throughout our simulations. Finally, we use a finite difference method to solve the time-dependent 
equations (1-4) in the QW plane. 
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Figure 1. Ambipolar diffusion coefficients vs. plasma density for an 8 nanometer InGaAs quantum well. 

First, we show the ambipolar diffusion coefficient for the InGaAs quantum well as a function of plasma density 
and temperature in Fig. 1. Here we only illustrate the nonlinearity in the diffusion coefficient and focus on its 
understanding from the physical point of view. Overall, the coefficient increases with plasma temperature, but in 
the intermediate density range many-body effects reduce it by as much as 20% in magnitude. At low density where 




classical statistics applies, both the plasma thermal energy term and the carrier-LO phonon scattering rates become 
density-independent and linearly proportional to the plasma temperature, whereas the bandgap renormalization term 
is negligible. By contrast, when the carriers become statistically degenerate at high density, the thermal energy term 
becomes the dominant contribution and increases with density linearly, whereas the scattering rates weakly decrease 
with density due to screening effects. As a result, the diffusion coefficient becomes a nonlinear function in the plasma 
density around the transparency point ( 1.9xl0 12 cm*" 2 in this study). Later in this section, we will show that this 
nonlinear diffusion coefficient Dj\r N can be substituted with an equivalent constant one as far as DC operations and 
AC-modulation performance are concerned. 



Figure 2. DC effects of lateral diffusion: (a) Near-field laser intensity vs. relative injection level J/Jth — 1 for four 
cases. Case (F): solid curve with circle symbols for nonlinear D NN , and cases (a)-(c): constant D N N s; (b) threshold 
current (circle symbols only) and slope efficiency (solid curve with circle symbols) vs. diffusion coefficient. The solid 
curve through the threshold current data points is a linear fit to the data. The equivalent constant coefficient data 
have been added (at 34.34 cm 2 /s). Refer to Figure 5(a) and text for further explanation. Denoted with a plus symbol 
is the nonlinear diffusion result. 


Next, we present the numerical DC operating characteristics of the VCSEL. The results are presented in a scaling 
manner where the relative injection level is adopted as the scaling parameter. As mentioned, the field intensity is for 
the center point of the device so that a constant conversion factor is required to obtain the far field intensity, and 
so does the slope efficiency. Nevertheless, this conversion does not affect our conclusions. Throughout this work, we 
choose four cases with representative diffusion coefficients to show the effects of lateral diffusion. Case (F) corresponds 
to employment of the full nonlinear diffusion coefficient Djvn(N) in simulation, case (a) features a constant coefficient 
of 54.115 cm 2 /s wffiich corresponds to a plasma density of 2.5xl0 12 cm -2 in relation with D nn (N), case (b) has 
23.838 cm 2 /s which corresponds to 7xlO n cm' 2 , and case (c) has 7.5 cm 2 /s which corresponds to 2.5 xlO 10 cm -2 . 
We opt to work with these coefficients because (1) the nonlinear case (F) falls among them, for both DC and AC- 
modulation characteristics; and (2) the range of values for the coefficients covers typical corresponding carrier denshy 
where appreciable change in the nonlinear coefficient is observed in Fig. 1. Therefore, the effects of lateral carrier 
density diffusion are demonstrated through comparison of the four cases. In Fig. 2, we plot the near field intensity 
as a function of the relative injection level J/Jth — 1 in Panel (a). As seen, both linear and nonlinear cases scale 
linearly with the relative injection level. This scaling behavior is referenced when we discuss the AC effects of lateral 
diffusion later. To better illustrate the effects, the threshold current (circle symbols only) and slope efficiency (solid 





curve with circle symbols) are plotted as a function of diffusion coefficient in Panel (b). Simultaneously shown m the 
plot are a linear fit to the threshold current and data points of the equivalent constant coefficient case (with a value 
of 34 34 cm 2 /s) for the nonlinear case (F). We shall discuss in detail the nonlinear case and its equivalent linear case 
later when the related result is presented in Fig. 5. Therefore, it is demonstrated that the induced increase in the 
threshold current by lateral carrier diffusion is linearly proportional to the diffusion coefficient. Correspondingly, a 
decrease in the slope efficiency results. 
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Figure 3. (a) AC-modulation spectral response at constant relative injection level of 20% for the same diffusion cases 
as in Figure 2; and (b) AC-modulation spectral response at constant difference injection level of A J - 0.312 K A/cm 
for comparison. 


As first seen from the DC effects of lateral diffusion presented above, the relative injection level seems to be a 
good controlling parameter for identifying the static and dynamic nature of the AC effects of lateral diffusion. Shown 
in Panel (a) of Fig. 3 are numerical results for the spectral response of the VCSEL under AC-modulation at the 
same relative injection level of 20% for the four diffusion cases. For comparison, in Panel (b) we plot the results at 
the same difference injection level of A J = J - Jth = 0.312 KA/cm 2 . It is obvious that to some extent the results 
comparatively exhibit the expected scaling behavior for the spectral response when the AC effects of lateral diffusion 
is purely of static nature. However, due to the associated change in the slope efficiency with the threshold current, 
the evidence is not overwhelmingly convergent upon a clear identification. Therefore, we take these results as a 
strong hint for the purely static nature of the AC effects of lateral diffusion and move on to performing a hneshape 
analysis on the spectral response of the VCSEL. 

Despite the fact that there is no known analytical guidance for lineshape analysis on the spectral response of a 
VCSEL under AC-modulation with diffusion process included, the suggestion comes from two sides. First, mrea l y, 
it is a common practice to model the spectral response using the classical linear rate-equation approach. Doing so 
reveals that the diffusive process affects the device physics of the VCSEL in such a manner that both the scaling 
behavior and spectral response are somewhat intact in functional form. Furthermore, the scaling behavior presented 
above provides additional support for a lineshape analysis since different responses tend to collapse together. The 
essence of a lineshape analysis lies in its ability to identify resonance position and dephasing mechanisms. In this 
work, we do not attempt to evaluate the latter. Instead, the analysis is used to find the former together with the 
peak amplitude and the results are checked to see whether they scale with the relative injection level. We argue that 
their scaling behavior would correlate with the DC effects of lateral diffusion and thus provide a strong but indirect 
evidence for identification of the static nature of the AC effects of lateral diffusion. 


Now, we are ready to present such evidence. In Panel (a) of Fig. 4 we show spectral response results for the evice 
at different injection levels but at a constant diffusion coefficient of 23.838 cm 2 /s. One fit curve (dot-dashed curve) 
to the case with the largest injection level (dashed curve) is shown. The other fits similarly show agreement near 
resonance. The results tell a few things: First and foremost, the functional form of the lineshape is indeed intact 
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Figure 4. (a) AC-moduIation spectral response at different injection levels and under a constant diffusion coefficient 
of 23.838 cm 2 /s. The dot-dashed curve is a lineshape fit to the long dashed curve, as explained in the text; and (b) 
Resonance frequency (plus symbols only) and resonance amplitude (circle symbols only) vs. relative injection level. 
The solid curves are linear fit to the data points, respectively. 


and given by 1/a; 2 — cu# + ia; 7 . 19 Second, diffusive process does cause growing deviations from this standard form at 
both low and high frequency sides. The nature of these deviations can be either static or dynamic or both, which is 
a subject of future study. Last, we point out that without an adequate assessment of the effects of diffusion, analysis 
performed based on the classical linear rate-equation model is expected to bear systematic discrepancy. Next, we 
show the resonance position and peak amplitude obtained from lineshape fitting as a function of the square root of 
the relative injection level in Panel (b). Note that in the plot the solid curves are linear fit to the data points. As seen, 
the resonance frequency scales with y/JJ Jth ~ 1 and the peak amplitude linearly depends on the relative injection 
level. The same scaling behavior is observed for the other diffusion cases, which are not shown. The scaling behavior 
of the resonance frequenc} 7 is expected for a coupled two-component system where one component experiences gain 
from excess above threshold of the other component. Therefore, results shown here present the evidence we need to 
reach our conclusion that the AC effects of lateral diffusion is purely of static nature, which is comprehensible from 
the physical point of view that diffusive process occurs at a slower timescale than the about 10 GHz modulation rate. 

It is interesting to take a closer examination of the classical linear rate-equation model and come to a more 
penetrating understanding of the effects of lateral diffusion. Basically the model 1 ^ features a linear gain approximation 
and no diffusion. For better discussion, we list the spectral response solution of the model below; 


p(u) = 


v 2 - iv( 7 at + AVo) - kAVo ’ 


( 6 ) 


where 8J is the amplitude of AC-modulation current density, A is the linear differential gain, and Vo is the DC laser 
fi eld intensit y. The resonance frequency is given from the above equation as u R = -JkAV 0 = y /'y N (n + ArN tr ) x 
V J/Jth - 1 with N tr being the transparency density of the VCSEL. There are a few ways that lateral diffusion 
alters this description. As revealed in Panel (a) of Fig. 4, deviations are present at both low and high frequency 
sides. Furthermore, as demonstrated, the scaling behaviors of the resonance frequency and peak amplitude are still 
valid, but their scaling prefactors are diffusion-dependent, which shall be commented upon in the next concluding 





section. Also, analysis based on 
the effects of lateral diffusion. 


the classical model surely introduces systematic errors without taking into account 






Cl 

CD 

TD 

C 

o 

"•4— * 

CO 


=3 

T3 

o 


£ 



5 15 25 35 

U ( GHz > 


Figure 5. Comparison of (a) near-field laser intensity vs. injection level and (b) AC-modulation spectral response at 
the same relative injection level of 20% for the nonlinear diffusion coefficient case [(F): solid curve with plus symbols] 
and its equivalent constant diffusion coefficient case [(const. D NN ): long dashed curve with circle symbols]. 


Last we show that an equivalent constant diffusion coefficient exists to reproduce the nonlinear diffusion char- 
acteristics by presenting a comparison of the DC and AC-modulation performance between the nonlinear diffusion 
case (F) with such a constant coefficient case. The coefficient for the latter case is chosen as follows: We first obtain 
the L-I curve for case (F) and find the corresponding threshold injection current. Then using the linear fit to the 
constant diffusion coefficient cases, as shown in Panel (b) in Fig. 2, we find the corresponding constant diffusion 
coefficient for the threshold current in case (F). This process is indicated in the figure by dashed lines with arrows 
pointing out the procedural flow. The simulation results are shown comparatively in Fig. 5. Needless to say, the 
equivalence is established. There are a few points worth mentioning: First, we have given quantitative evidence 
for the equivalence and an unambiguous criterion to locate the value of the equivalent constant coefficient Sec- 
ond, this criterion is quantitatively different from previous studies 9 - 11 which pointed to an average constant value 
of the diffusion coefficient’. In this work, the average value, if given at half the pinned plasma ensi y ( pin 
1.893 xlO 12 cm- 2 ) in the middle of the VCSEL, is 29.22 cm 2 /s, which is about 15% below the equivalent coefficient 
of 34.34 cm 2 Is. Certainly, other averaging procedures could be introduced to reconcile with this offset, but our crite- 
rion is physically meaningful and straightforward. However, the reality is that the linear fit to the constant diffusion 
coefficient cases is nonexistent. Therefore, it brings out the last and most important point we would like to make, 
i.e the nonlinearity in the diffusion coefficient is important in a self-consistent manner Seemingly the equivalence 
of the nonlinear diffusion coefficient to a constant one should deny the necessity of a full nonlinear description of 
the diffusive process. We emphasize that this conclusion is restricted to the studied range and can not extrapolate 
without further investigation. Ramifications could result in other complex scenarios. Even in the restricted range, 
this equivalence has been established numerically through comparison. Therefore, it is impossible to iden l y e 
equivalent coefficient without first characterizing the nonlinear case. 





4. SUMMARY AND CONCLUDING REMARKS 

In summary, we have presented the microscopically derived density diffusion coefficient as a function of plasma 
density in the AD regime. The effects of lateral diffusion on DC operation and AC-modulation characteristics are 
studied through comparison of the nonlinear diffusion coefficient case with three constant coefficient cases and its 
equivalent constant coefficient case. Specifically, a scaling scheme and lineshape analysis are deployed in identification 
of the static and dynamic nature of the spectral response of the VCSEL under AC-modulation. The last numerical 
result presented is a direct comparison of the DC and AC-modulation performance between the nonlinear diffusion 
case and its equivalent constant diffusion coefficient case. 

Before closure of this paper, a few remarks are given. First, the nonlinearity in the diffusion coefficient does alter 
performance of the device. Second, despite the fact that we have shown that such a nonlinearity can be modeled 
with an equivalent constant coefficient, the answer to how to find the equivalent value is still pending. We hint 
here on how to proceed. As shown, the increase in the threshold current is proportional to the coefficient. Given 
that the proportionality coefficient is known, empirically or analytically, the equivalent constant coefficient is thus 
determined. By analyzing the AC-modulation response, if some analytical guidance were supplied, in principle we 
then should be able to obtain such critical device parameters as the differential gain, carrier lifetime, and photon 
lifetime. So we see a critical need for development of such an analytical understanding of the device with the effects 
of lateral diffusion included. 

In conclusion, we have shown that within the studied range, i.e., 50% within injection current threshold for 
a single mode gain-guided VCSEL, not only the nonlinear diffusion coefficient can be replaced with an equivalent 
constant, but also a numerical scheme is provided to locate the constant. The numerical criterion established is that 
both coefficients yield the same threshold current. For the DC effects of lateral carrier density diffusion, it is found 
that the increase in threshold current depends linearly on the coefficient, whereas the slope efficiency decreases with 
diffusion coefficient. For the AC effects of lateral carrier density diffusion, it is found that, the spectral response 
follows the lineshape function of 1/uj 2 - u>% + iu 7 near resonance and deviates at low and high frequency sides. Both 
the resonance frequency uir and resonance amplitude scale with the relative injection level J / Jth — 1* Furthermore, 
using the scaling scheme and lineshape analysis, it is concluded that the AC effects of lateral carrier density diffusion 
is purely of static nature and are reflections of the corresponding change in the threshold current and slope efficiency. 
Finally, we emphasize the importance of taking into account nonlinearity in the diffusion coefficient at a self-consistent 
level. 
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The lineshape is determined by various interactions, such as that between laser field and plasma density, w ic 
gives the resonance behavior, and those incoherent ones between the coupled field-plasma system and t e 
ambient, which produce the linewidth. In this work, vertical carrier transport and capture is not considered so 
their corresponding contributions to the spectral response is missing, which is conventionally modeled with a 
multiplicative factor 1/1 +iur tc where r tc is a lumped time constant accounting for the corresponding relaxationa 
mechanisms. 




